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1. Introduction 

Various methods are used to solve systems of ordinary differential equations (SODEs) in 
mathematics and applied sciences. The most popular ones include factorization |15j . linearization 
perturbation [141120] . closure approximation mm, discretization [23j and Adomian methods 
H H H [H El] to cite a few. Recently, in an interesting paper [3, Bagarello developed a non- 
commutative method, based on the quantum mechanics formalism, for the analysis of systems of 
ODEs. He provided solutions of some systems described by unbounded self-adjoint and densely 
defined Hamiltonian operators and discussed corresponding integrals of motion. 

The possibility of acting, in an appropriate way, on systems of ODEs governing physical 
phenomena or other mechanisms in nature gives rise to the control theory. Optimal control 
problems can be solved by the so-called Pontryagin’s method. 

One of the most important systems encountered in the literature is certainly the system of 
harmonic oscillators, widely used as a basic tool in physics Esuni- This is motivated by their 
role in many applications in various fields of physics and technology. The harmonic oscillators 
describe a wide class of physical models. Besides, their properties are well known. This explains 
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why they are often used as a preliminary tool in order to gain initial insight into complex systems. 
For instance, quantum mechanics as well as optimal control were first illustrated using systems of 
harmonic oscillators. A typical controlled harmonic oscillator is a pendulum made of a string and 
a ball moving in a vertical plane and subjected to a force (also called the control) whose role is to 
stress the system to rest in a minimum time. 

Recently, Andresen et al [4] invesigated the control of an oscillator using the frequency vari¬ 
ation. Dasanayake m and Van Dooren [25] . (see also references therein), performed numerical 
solutions of optimal control problems by pseudospectral methods and Chebyshev series. 

In the present work, for the first time to the best of our knowledge, the Bagarello’s noncom- 
mutative approach is applied to the control of systems of nonlinear ordinary differential equations. 

The paper is organized as follows. In Section 2, we present the considered optimal control 
problem. We derive the state and adjoint systems of ordinary differential equations using the 
Pontryagin’s maximum principle [T^[2T]. Then, we use the Bagarello’s operator method to solve 
such systems. In Section 3, harmonic oscillator systems are controled and discussed. 


2. Theoretical framework 


We consider the following optimal control problem mm- 


x = f{x,u), x{0) = Xo 

x{T) = xi 
uGV C 

r 

min / f°{x,u)dt, 

“eK 7o 


( 1 ) 

( 2 ) 

(3) 

(4) 


where x,Xo,Xi € A"; D V 3 u is the control, V the control set; T stands for the moment 
and E is the Euclidean space. All the functions occurring in the formulation of the problem are 
assumed to be differentiable with respect to {x,u). Further, the minimum should be in a class of 
measurable and bounded functions u, defined on [0, T]. For instance, V can be represented in the 
form 


V = {u G E’’ : hi{u) < 0, f G Ii;hi{u) = 0, f G h} 


where the functions hi express the constraints on the control; Ii and I 2 are finite disjoint index 
sets. 

Let us introduce the Hamilton-Pontryagin’s function 


r]{x, u, Ip) = ipf*{x, u) - tp^fix, u) 


where / = (/i,..., /„) G E”; E” 3 tp = {ipi, ..., tpn) is the adjoint variable depending on t; * 
stands for the transposition operation. ■0° is a positive number. 

By virtue of the Pontryagin’s maximum principle [mill], the existence of adjoint functions 
is a necessary condition to link the control system of ordinary differential equations to a functional 
objective. The optimal control u* is then characterized in terms of the state and adjoint functions. 

Thus, given an optimal control u* and the corresponding system o, there exists an adjoint 
variable tp satisfying the following equation: 


dip 

dt 


There results the Cauchy problem 


dr], , 


{x,ip),ip). 


X = f(x,U*(x,'lp)), x(0)=xo 

Ip = -^{x,U*{x,1p),1p), 0(0) =00 


with the unknowns x,ip. 


Setting y 



/ f{x,u*{x,ip)) 

\-^{x,u*{x,ip),ip) 



(5) 

( 6 ) 

(7) 
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we can write the latter in the form: 


or, explicitly, 


y = 

= 9{y), 2/(0) = 


2/1 

= 51(2/1,2/2,- 

■-,52n) 

2/2 

= 52(51,52,- 

■-,52n) 

2/2n 

= 52n(5l,52, 

---,52n) 


5i(0) = 5°, 

, 5 = 1,- 


( 8 ) 


(9) 


We suppose that the functions gj are such that the solution of the Cauchy problem (l9|) exists and 
is unique. 

In this work, we intend to solve this system by Bagarello’s approach ([5, and references 
therein), based on a quantum mechanical formalism. Given a quantum mechanical system S and 
the related set of observables Os, i- e., the set of all self-adjoint bounded (or more often un¬ 
bounded) operators describing S, the evolution of any observable Y S Os satisfies the Heisenberg 
equation of motion (HOEM)[5: 

jY{t)=z[H,Yit)]. (10) 

Here [^, H] = AB — BA is the commutator between A, B G Os',H is assumed to be a densely 
defined self-adjoint Hamiltonian operator of the system acting on some Hilbert space given by 


m 


2n 


^i9o) = 9 51 {Pidjivi, vD + 9jiyi,y^, ■■■, yDpj} 


( 11 ) 


i=i 


where go = (^gi{Y°), g 2 (Y°),..., g 2 n(Y^)^ ,Y° = (yj, ..., y^„); the initial position y° is con¬ 
sidered as an operator acting on the Hilbert space 'H, and pj is its canonical conjugate momentum 
operator such that 


[y°,Pk] = iSj,kl, j,k = l,...,2n 

[y^yk] = [Pj,Pk]=0, j,k = l,...,2n. 

Standard results in quantum mechanics show that for any differentiable functions 
P {Pi, yl,---, vln) > and p {P 1 ,P 2 ,. •. ,P 2 n), we have 

[pj,v>iyi,y 2 ,---,yL)] = j = i,...,2n 


( 12 ) 

(13) 


[yj,f’{Pl,P2,---,P2n)] = j = l,..., 2 n. 


(14) 

(15) 


Theorem 1. [T] // the functions gj are holomorphic, a formal solution of the HOEM (0) is 

r(<) (16) 

where Y^ is the initial value of Y (t) and H does not depend explicitly on time. 

Furthermore if H is bounded, we get 

iitf 




fc >0 


k\ 


-[H,Y^ 


(17) 


where [A,B]k is the multiple commutator recursively defined as : 

[A,B]o=B; [A,B]k = [A,[A,B]k-i]. 

The concepts of integral of motion and extended integral of motion of the system 
defined as in [ 7 ]. 


are 
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Definition 1. Any holomorphic function I depending on the variables yj, such that 

.. ,y 2 n{t)) = Iq Vi, Jo = const., 
is called an Integral of motion (loM) of system 0). 

Definition 2. We call extended integral of motion (EIoM) of the system any holomorphic 
function J depending on the yj,Pj such that 

Jiyi{t),y2{t),...,y2n(t),pi{t),P2{t),...,P2n(t)) = Jq Vt, Jq = COnst., 
where yjit) = pj{t) = e^^*pje~'^^^, j = 1 , ..., 2 n. 

One of the main advantages of the Bagarello’s strategy is to provide at hand a good approxi¬ 
mation of the solution of the SODE. The th order approximation (iV G N) is 

YN{t) = Y.^[H,Y%. (18) 

k=0 

Using the integral of motion, we can estimate the approximation’s error of the unknown exact 
solution of the derived Cauchy problem. In this case, one can proceed as follows : let Y]sf{t) be 
the approximated solution of Y{t). We may thus compute and evaluate the error 

AN{t)=I{YM{t))-I(Y°). (19) 

This approach is useful for our study because formula (fTSl) provides a simple approximation 
scheme for exact solutions to ordinary differential equations. Further the estimation of the error 
Ajv(t) tells us how good this approximation is. 

In the next section, we apply this approach to time-optimal control problems. 


3. Applications 


3.1. Problem 1. In classical mechanics, we study the forced harmonic oscillator. The cor¬ 
responding optimal control problem has the form [18] 


Xi 

X2 

a;i(0) 

Xi{T) 

-1 

T 


= X2 

= + u{t) 

= a:?; X2{Q)=xI-, 

= 0 ; X2{T) = 0 ; 

< u<l 
—>inf. 


( 20 ) 


Here t and T denote the time, x = {xi,X 2 ) € E'^, u is the control and uj the oscillation frequency. 

In this problem, the infimum is sought in a class of controls u(t),t > 0. We study a control u 
that moves the point a:(0) to the point x{T) — 0 in accordance with the corresponding solution of 
(I20p during the time T. Then we solve the boundary value problem derived from the Pontryagin’s 
maximum principle. 

Let us consider the Hamilton-Pontryagin’s function 


ri{x, U, Ip) = -1 + 1plX2 + ip 2 {-‘x'^xi + u) 


( 21 ) 


where ip = {ipi,ip 2 ) G E^- 

According to Pontryagin’s maximum principle, we obtain the control u in the form 


f 1 if '02 (t) > 0 

( -1 if -02 (t) < 0. 


Now we solve the Cauchy problem derived from the Pontryagin’s maximum principle using the 
technique developed by Bagarello. 
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In the case u{t) = 1, we can write 


Xi = X2 

X2 = —Uj'^Xi + 1 
Ipl = Uj'^1p2 
'Ip2 = -l/’l 

xi(0)=a;?; 0 : 2 ( 0 ) =^ 2 ; V'i(0) = V'i; V' 2 ( 0 )=V' 2 - 

Introducing the new variable ii = xi- Ti rewrite (l2^ in the form 


( 22 ) 


Xi = X 2 
±2 = —Up'Xi 
Ipl = 0j‘^'tp2 
'02 = -i’l 


(23) 


In the case u{t) = — I we get 


or, equivalently, denoting xi = xi 




Xi = 

X2 

±2 = 

—u'^Xi — 1 

Ipl = 

01^02 

'02 = 

-01 

( Xi 

= X2 

J X2 

= —Uj'^Xi 

] 

= w ^02 

[ 02 

1 

II 


(24) 


(25) 


with the previous initial conditions . 

Equations ([23|) and ([25|l reduce to the equivalent system 


yi = y2 

m = 


^3 = uj'^yi 


(26) 


2/4 = -2/3 

where yi = xi or y^= xi,y 2 = 0 : 2 , 2/3 = '01,2/4 = 02- 

Here, instead of using standard methods widely spread in the literature, we show that Bagarello’s 
suggestion can also be exploited to solve this optimal control problem and to reproduce known 
results. Indeed, using the definition of the commutator 

9]P3 = - [Pj : 93 ] +Pj93, J = 1, ■ • ■, 4, 
we compute the Hamiltonian of the previous system in the form 

H = Piyl - uj‘^P2yi + w^P32/4 - PiPs- 
Developing the solution of (E51) as an infinite series: 


(27) 


fyi{t)\ 


fy°\ 

2/2(0 


2/2 

2/3(0 


2/3 

\yiit)J 


w 


it 




(yT\ 

2/2 
2/3 

\yV\ 


2 ! 


H, 


fyi\ 

2/2 

2/3° 


3 ! 

H, 

fy°i\ 

2/2 

2/3° 


\yV 


2 


\yV 



we calculate the multiple commutator 


f[H:y°i]2A 

[H.y°2]2, 


u3^Py° 


[H.yl]2, 


LO^Py° 

7 

\[H.yi]2j 


W^y2) 



([H.y'i], 

[H,yV\, 

[H,yV\, 

V[^,2/4]: 


2p+l^ 


/ -iuj^Pyl \ 

2p+l 



2p+l 

2p+l/ 


—iLtP‘P^'^y\ 

\ / 


p = 0,l,2,3,... 
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and obtain 

/ \ 

Ui cos(a;t) H —- sin(w<) 

(jJ 

j /2 cos(wt) — UiUJ sin(a;t) 
j/g cos{ujt) + y^uj sin(ajt) 

y cos{ujt) - ^ sin(tLit) / 

In order to determine the time T such that 

/ xi(r) = o 
1 X2iT)=0 ’ 


y2{t) 

ysit) 

\ 2 / 4 (t)/ 


and taking into account yi = xi = xi -y 2 = 3^2, we obtain 

UJ^ 


Xi{t)\ _ ({x°- cos(a;t) + ^ sin(a;t) + 


I - I ' UJ ' ' ' \JJ ' ' \JJ 

X2\t)) y X2 cos{ujt) — {xi — ■^)uj sm{ujt) 


and 


Thus 




tanwT = 


^{xi - -4) 


or ujT = arctan 


X 2 

~ ZJ?) 


T = — arctan ■ 

UJ 


^0 

X 2 


^{xl--^) 




Analogously, for yi = xi = Xi -\ -y 2 = 2 : 2 , we get 


= ({^1 + 1^) cos(a;t) + ^ sin(a;t) - ^ 


and the relations 


. . I - I \ UJ / ' ' UJ ' ' u 

X2\t)J \ 2:2 cos{ujt) — (a:° + a;sin(a;t) 


cos(wT) =H-^)sin(wT); sin(wT) = 


X 2 /UJ 


ojHxU:^)^ + ix°2)^ 


implying 


tanwT = 


X 2 




or ujT = arctan 


+ ZJ2-) 


f = t arctan 2:2 > 0- 

a; io[xl + ^y n 

Therefore, the optimal time solution to problem (l20l) is given by considering the following two 
cases : if 2:2 > 0, then T is the optimal time and if x® < 0, then T is the optimal time. 

In the sequel, we deal with more complicated problems whose solutions can be obtained by 
perturbative approaches. In this case, we show that as a possible candidate, Bagarello’s formalism 
can be easily implemented in a suitable and solvable form. 
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3.2. Problem 2. Let us consider the optimal control problem of a pendulum with large 
oscillations in the form 


Xl 

±2 

a;i( 0 ) 

MT) 

-1 

T 


= X2 


= —(3x2 — sina:i + u{t) 


= a;2(0)=a:2' 

= 0; X2{T) = 0; 

< R < 1 

—>■ inf 


(28) 


1 


u(t) = — F{t), 0 < t < r, is the control subjected to the following constraint : 
m 

u G 1^ = {m £ : |u| < l} ; 

Xi,X 2 are given positive constants, m is the mass, F the force and /3 > 0 a constant. 
The Hamilton-Pontryagin’s function is written as 

'q{x, u, V') = —1 + 'ipiX 2 + ip2{—!3x2 — sincci + u) 

where ij) = {tjji,tjj 2 ) G 

By virtue of the Pontryagin’s maximum principle, the control is given by 

1 if ip2it) > 0 

-1 if i>2{t) < 0. 

The associated Cauchy problem reads 


u{t) = 


Xi = X2 

X2 = — 13 x 2 — sinxi + u{t) 

-Ipl = 1p2 COS Xl 

1p2 = -V'l + /3V'2 


with the initial conditions 


a;i(0) = > 0; 0:2(0) = cca,> 0; V’l(0) = V'l; V'2(0) = V'2- 

Equivalently, this can be re-expressed as follows: 

yi = 2/2 

3)2 = -I 3 y 2 - sinj/i -I- u(t) 

2/3 = 2/4 cos yi 
m = -2/3 + Pyi 

2/i(0)=y?; 2/2(0) =2/°; 2/3(0) =2/g; 2/4(0) = 2/2 

where i/i = xi, 2/2 = 0:2,2/3 = 2/’i7 2/4 = V'2- 

The Hamiltonian of the system takes the form: 

H =Pi2/2 FP2 [-I 3 yl - sin2/? + u{t)] -h P32/4 cos 2/? +P4(-2/3 +/32/4)- 
The solution (fT^ becomes in this case: 


/2/?\ 


Y2{t) = 

where the commutators are given by 


2/2 

2/3° 

\yV 


+ it 


H, 


(yl\ 

2/2° 

2/3° 

\vV\ 


2 ! 




(y"\ 

2 / 2 ° 

2 / 3 ° 

\yV 


(29) 


(30) 


(31) 


(32) 


(33) 


/ 

H,y°i 

\ 


( -* 2/2 

\ 


H,y°2 


-i [-Pyl - sin 2 /? -f 

w(t)] 


'H.yt 



-iyl cos 2/1 

V 

H.yl 

} 


\ -*(-2/3 +/32/4) 

/ 
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/ 

H,y°' 

2 ] 

/ 


H,y°2 

2 



H,y°3 

2 


V 

H,y2 

2 ) 

V 


Finally we get 

• For u{t) = —1, 


yi.2{t) 

1/2.2 (t) 

and 

X2.2{t) 


{-if [-/32/2 - sinyi + u{t)] 

{-if {-/3 [-( 3 y 2 - sin7/5* + u{t)] - 008(7/5*)} 

{-if [- 2 / 2 I /4 sin(7/5') + (-7/g + I3yf) cos(7/5’)] 
i-if [- 2/4 cos(7/5’) + /3(-7/§ + Py^)] 


yi + tyl + ^{-Pyl - sin 7/5* - 1) 




v2/2 + - siny)* “ 1) + V [-/^(“/^l/” “ siny)* - f - vl cosy)* 


Xi + tx 2 + ^{—13x2 — sina;? — 1) 

C 2 + t{—j3x2 — sinX]* — 1) + ^ [—/3(—/3a;2 — sinx)* — 1) — a;® cos^J] 


For u{t) = 1, 


yi.2{t) 

y2.2{t) 

and 

Xl.2{t) 

*2.2 (t) 


y)* + tyl + ^{-Pyl - siny)* + 1) 

v 2/2 + t{-l3y2 - siny)* + 1) + ^ [-I3{-I3y2 - siny5’ + 1) - y 2 cosy)*] 


Xi + tx 2 + ^(—/3a;2 — sina;)* + 1) 

.3^2 + t{—Px 2 — sin a:? + 1) + ^ [—/3(—/Jx® — sinx? + 1) — x® cosx)*] 


Taking into account the relations 


*1.2(T) =0 
*2.2(T) =0 


and 


*1.2(T) =0 
* 2.2 ( T ) =0 


we obtain 


(i) f = 


{— 13 x 2 — sinxj — l)(a;2 + Pxf) + xjx^ cosx° 

{—{ 3 x 2 — sinx° — l)(sinx° + 1) — (2^2)^ cosx? 

and 02 [(aiC2 — a2Cif + (a2&i — ai62)(&iC2 — &2C1)] = 0 


(ii) T = 


with Oi, &i, Cl and 02 , & 2 j C 2 being respectively the coefficients of the polynomials * 1 , 2 (i) 
and * 2.2 (i)- 

(—^X 2 — sinx)* + l)(a ;2 + I3x\) + xjx^ cosx)* 


{—(3x2 — sinxi + l)(sinxi — 1) — {x^f cosx^ 
and A 2 [( 2 I 1 C 2 — A 2 Cif + {A 2 B 1 — AiB 2 ){BiC 2 — B 2 Ci)'^ = 0 


with Ai,Bi,Ci and A 2 ,B 2 ,C 2 being respectively the coefficients of the polynomials 
* 1.2 (t) and * 2.2 (i)- 

{—(3x2 — sinxi — l)(x 2 + (3xi) + XiX® cosx? 

{—(3x2 — sinxj — l)(sinxi + 1) — (a;^)^ cosxj 

{—(3x2 — sinx)* + l)(x 2 + (3xi) + x^x^ cosx)* 

{—(3x2 — sinXi + l)(sinxi — 1) — {x^f cosxi 

then T is the solution of the problem. Otherwise T is the solution. 

Let us estimate the error using the following integral of motion 

/(xi, X 2 ) = X 2 +/3xi + / [sinxi(r) — 77(r)] dr; with/(x?, x®) = x® +/3x?. 

Jo 


Then 
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(i) for u{t) = — 1, 


A2(t) = /(ii.2,i2.2) - - 1) - —ccacosa;? 


A2(f) 


sin 


< — Xn + it 

- 2 ^ 


Xi + tx 2 + -^{—13x2 — sin a;? — 1) 


+ 1> dr; 


leading to 

(ii) for u{t) = 1, 




<1 

< R for t G 


0 , 


—4 + y^ 16 + 2a;2R 


, with K a sufficiently small positive number. 


A 2 {t) = I{Xi, 2 ,X 2 . 2 ) - I{Xi,X 2 ) =t{-smx° + l) - —X 2 COSX° 


/o 


sm 


+ tx 2 + -^{—ldx 2 — sin a;? + 1) 


-l>dT; 


AaW 


^2 

< -^^2 + 


yielding 




A2(0 

< K for t G 


0 , 


—4 + y^ 16 + 2a:2N 


, with K a sufficiently small positive number. 


3.3. Problem 3. Let us examine the optimal control problem when the state equation is the 
Van der Pol equation mm 

X + £x{x^ — 1) + a: = u{t) (34) 

such that the control u{t) G [a, /3] . The problem has the form 


Xi 

= X2 

X2 

= — Xi + 6X2(1 — x 0 + u(t) 

Xi(0) 

= x?; X2(0)=x^; 

Xi(T) 

= 0; X2(T) = 0; 

a 

< u < P 

T 

—> inf 


(35) 


where e, a, /3 are real constants. 

The Hamilton-Pontryagin’s function is 


■q{x, u, -0) = -1 + '4’iX2 + 02 [-a^i + £X 2{1 - xf) + u(t)] . 


(36) 


According to the Pontryagin’s maximum principle, the supremum of the function rj depending on 
xi, a: 2 , 01 , 02 , rt with respect to u is reached when the control takes the following form: 


u{t) = 


0 if 02(0 > 0 

a if 02 0) < 0 


Two cases are examined 

(i) Taking u{t) = a, we get the following Cauchy problem 


Xi = X2 

±2 = —Xi + £X2{1 — xl) + I 
01 = (1 + 2£XiX2)02 


02 = - [01 +e(l - a;i)02] 

^xi(0)=x?; X2(0)=x^; 0i(O) = 0?; 02(0) = 0°, 


( 37 ) 
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which, according to m, can be reduced to 

yi=y 2 

m = -yi + ey 2 {l -yf)+a 

< ya = (1 + 2 eyiy 2 )yi ( 38 ) 

y4 = - [y3 + e(l - y1)y4\ 

.2/i(0)=2/i; J/2(0) = 2/^; 1/3(0) = y§; 2/4(0) = 

The Hamiltonian of the system is computed as follows 

H = piyl + P 2 {-yl + eyl [l - (y?)^] + a} + ^ 3(1 + 2 eylyl)yl 

+Pi {2/3 + £ [1 - (1/?)^] 2/4} + *£ [1 - iyif] ■ ( 39 ) 

The first order approximation solutions are given by 


/ 

H,yf 

] 


/ -*2/2 \ 


H.y°' 



-i {-y? + ey^ [l - (y)’)^] + a} 


H,yl 



-i [(1 + 2ey)’yg)y2] 

V 

H,yi 



\ -*{2/3+£[l-(2/i)^]2/4} / 


{yi.i{t)\ ^ ( 2 /?+^ 2/2 ^ 

\y 2 . 1 it)) V 2/2 {- 2/1 +£ 2/2 [1 - ( 2 / 1 )^] +a}y ■ 

In the original variables 

_ f x'i+tx'^ \ 

\x 2 . 1 it)) \ai2 + i {“^1 + £^2 [1 “ ( 2 ^ 1 )^] + a}y ' 

(ii) For u{t) = {3, the corresponding derived Cauchy problem is put in the form: 


±1 = X 2 

±2 = -Xi + 6 X 2(1 - xf) + P 

< ipi = (1 + 2exiX2)ip2 (40) 

ip2 = - [tpi + e{l - xl)ip2] 

^xi(0)=x^; X2(0)=a:^; i/!i(0) = V'?; 1 ^ 2 ( 0 ) = ipt 


This SODE differs from (l38l) only by the term /3 replacing a. Then, replacing mu- 
tatis mutandis a by /3, the SODE remains the same as in (1381) . Hence, the first order 
approximation solutions are given in the form: 

yi+ty2 

\y2.i{t)) V2/2 +1 {-2/1 + £2/2 [1 - (2/1)^] + / 3 } 


or, equivalently, in terms of the original variables 
fxi.i(t)\_f x^+tx^ 

\X 2 . 1 it)) +t{-x°i + 6 x 1 [1 - ( 2 ^ 1 )^] +/3} 


Taking into account the relations 


rii,i(r) = o rxi,i(T) = o 

\ i 2 . i ( r ) = o \ x2.i(t)^q 


we get the following results 


T = 


x\ — 6 X 2 [1 — (ic?)^] — a 
with the relations 


> 0 


and 


T = 


,0 [1 _ (^0)2] _ p 


> 0 , 


( 2 ^ 2 )^ + ( 2 ^ 1 )^ ~ 6 X 1 X 2 [1 — ( 2 :?)^] — axi = 0 or 


(^0)2^ (^0)2 _ ^^0^0 
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The optimal time solution to the problem is 

f if 


0-2 


< 


xl 


min(T, T) = <^ “ £^2 (1 — — a x\ — £X 2 (1 — — P ’ 

T otherwise. 

An integral of motion of the SODE (l3^ is 

J(a:i,a:2) = X2 + £ (^^x^ - Xi^ + J [xi{t) - u{t)] dr 


J{xi,X2) = + £ { ^{x'i f - x'i 


Setting 

(i) for u{t) = a, 


1 


Ai(t) = J{xi.i,X2.i) - J{x°,X2) =-£ix2)^t^ + 


x'-^ + £xl{x^f 


and using the Cardan’s formulas for cubic polynomials, we obtain 
Ai(t) < K for t G 


+ V 2"^“ 3 


A 


with K a sufficiently small positive number and 

,3 /n\2 1 


« = (I) +(l) ■ ■'= 27 '"’-'®' 


A = 

(ii) and for u(t) = /3, 


1 + 3 \£XiX2 


B = 


3k 


e kEn 


£ Xn 


|3 ’ 


Ai(t) = -e(x“rt^ + 


+£ X ?( x ^)2 


we obtain 


Ai(t) 


< K for t G 


0, + V^+ ^ 


with K a sufficiently small positive number and 

,3 x2 1 


" (?) +(l) ’ ^'""3^1’ 

3 

Ai = 


2+3|£x?x0| 3^ 

- , ,, -1 — , ,, „ 


£ Xo 


£ Xn 


3 ■ 


(41) 


4. Concluding remarks 


In this paper, we have successfully extended the domain of applicability of Bagarello’s ap¬ 
proach, developed for the analysis of ordinary differential equations, to optimal control problems. 
The relevance of the derived approximation scheme (1181) and error estimation given in formula 
HU) has been exploited to investigate time-optimal problems of forced harmonic oscillator systems. 
Three particular cases have been explicitly treated and discussed. 
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